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Abstract 

We develop a formalism for particle production in a field theory cou- 
pled to a strong time-dependent external source. An example of such a 
theory is the Color Glass Condensate. We derive a formula, in terms of 
cut vacuum-vacuum Feynman graphs, for the probability of producing a 
given number of particles. This formula is valid to all orders in the cou- 
pling constant. The distribution of multiplicities is non-Poissonian, even 
in the classical approximation. We investigate an alternative method of 
calculating the mean multiplicity. At leading order, the average multiplic- 
ity can be expressed in terms of retarded solutions of classical equations of 
motion. We demonstrate that the average multiplicity at next-to-leading 
order can be formulated as an initial value problem by solving equations 
of motion for small fluctuation fields with retarded boundary conditions. 
The variance of the distribution can be calculated in a similar fashion. Our 
formalism therefore provides a framework to compute from first principles 
particle production in proton-nucleus and nucleus-nucleus collisions be- 
yond leading order in the coupling constant and to all orders in the source 
density. We also provide a transparent interpretation (in conventional 
field theory language) of the well known Abramovsky-Gribov-Kancheli 
(AGK) cancellations. Explicit connections are made between the frame- 
work for multi-particle production developed here and the framework of 
Reggeon field theory. 
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1 Introduction 



The study of multiparticle production in hadron collisions at high energies is an 
outstanding problem in the study of the strong interactions. With the advent 
of QCD, it was understood that semi-hard particle production in high energy 
hadronic interactions is dominated by interactions between partons having a 
small fraction x of the longitudinal momentum of the incoming nucleons. In the 
Regge limit of small x and fixed momentum transfer squared - corresponding 
to very large center of mass energies squared s — the Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) evolution equation |1I2| predicts that parton densities grow 
very rapidly with decreasing x. Because this rapid growth in the Regge limit 
corresponds to very large phase space densities of partons in hadronic wave- 
functions, it was proposed that saturation effects may play an important role in 
hadronic collisions at very high energies i3l4i5i6i . Saturation effects slow down 
the growth of parton densities relative to that of BFKL evolution and may 
provide the mechanism for the unitarization of cross-sections at high energies. 

The large phase space density suggests that the small x partons can be 
described by a classical color field rather than as particles |7I8I9| . More precisely, 
the McLerran-Venugopalan (MV) model proposes a dual description, whereby 
small X partons are described by a classical field and the large x partons act 
as color sources for the classical field. The original model considered a large 
nucleus containing a large number of large x partons (at least 3^ valence quarks 
where A is the atomic number of the nucleus). In this limit, they produce a 
strong color source and one has to solve the full classical Yang-Mills equations to 
find the classical field. This procedure properly incorporates the recombination 
interactions that are responsible for gluon saturation. In the MV model, the 
large x color sources are described by a Gaussian statistical distribution |7I1U) . 
A more general form of this statistical distribution, for SU{Nc) gauge theories, 
valid for large A and moderate x, is given in Ref. |11I12| . 

The separation between large x and small x, inherent to the dual description 
of the MV model, is somewhat arbitrary and has been exploited to derive a 
renormalization group (RG) equation, the JIMWLK equation |13I14I15I1HI17I 
118119120] . This functional RG equation describes the change in the statistical 
distribution of color sources with x. It can be expressed as an infinite hierarchy 
of evolution equations for correlators j21| . and has a useful (and tremendously 
simpler) large Nc and large A mean-field approximation |22| . known as the 
Balitsky-Kovchegov equation. This general framework is often referred to as 
the Golor Glass Condensate (CGG) 23 24 25j. 

The discussion above refers to the properties of hadronic wave functions at 
high energies. To compute particle production in the GGC framework, (a) one 
must find the distribution of color sources in the projectiles, either by solving 
the JIMWLK equation or by using a model such as the MV model, and (b) one 
must calculate the production of the relevant particles from the sources. In this 
paper, we will assume that the source distributions are known and shall focus 
on particle production in the presence of these sources. 

At leading order, one knows how to compute the production of quarks and 
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gluons from strong color sources. Inclusive gluon production is obtained by solv- 
ing the classical Yang-Mills equations for two color sources moving at the speed 
of light in opposite directions |26l27l28j . This problem has been solved numeri- 
cally in |29I3UI31I32I33I34] for the boost-invariant case. A first computation for 
the boost non-invariant case has also been performed recently The multi- 
plicity of quark-pairs is computed from the quark propagator in the background 
field of [2913013 Il32i^5l51| - it has been studied numerically in |36l37j . 

These results summarize the state of the art in studies of particle produc- 
tion from strong classical color sources. In the present paper, we will system- 
atically examine particle production in a field theory coupled to an external 
time-dependent classical source. In particular, we will consider how one com- 
putes results beyond the above mentioned leading order results. To avoid the 
technical complications of gauge theories, we shall focus our discussion on the 
theory of a real scalar field. The complexities of QCD, such as gauge invari- 
ance, are of course extremely important. However, several of the lessons gained 
from the simpler theory considered here should also apply to studies of particle 
production in QCD. These will be applied in follow ups to this paper. 

This paper is organized as follows. We begin, in the next section, by setting 
up the model. We recall well known results for counting powers of the coupling 
constant and of the powers of h; in particular we discuss how these are modified 
for field theories with strong external sources. We discuss the vacuum- vacuum 
diagrams that will play a crucial role in the rest of the paper. We show clearly 
that the sum of the vacuum-vacuum diagrams is not a trivial phase (as is the 
case for field theories in the vacuum) but gives a non-trivial contribution to the 
probability in theories where the field is coupled to a time-dependent source. 

In section 13 we use Cutkosky's cutting rules to derive a formula for the 
probability P„ of producing n particles, in terms of cut vacuum-vacuum dia- 
grams. This formula, for theories with strong sources, is new and is valid to 
all orders in the coupling constant. While the formula at this stage is not use- 
ful for explicit computations, it shows unambiguously that the distribution of 
the multiplicities of produced particles is not Poissonian, even in the classical 
approximation. We also derive a compact formula for the generating function 
F{x) of the probabilities P„. 

In section 0] we discuss in detail the calculation of the first moment of the 
distribution, namely, the average multiplicity (n) = X^n'^-^"- We obtain a 
formula which is valid to all orders for ^n). We work it out explicitly at leading 
and next-to-leading order. It is well known that the leading order expression 
can be expressed in terms of classical fields that are computed by solving partial 
differential equations, with retarded boundary conditions, for the classical fields. 
What is remarkable and a major result of this work is the fact that one can 
similarly compute the multiplicity at NLO from the small fluctuation flelds (in 
the presence of the classical background field) with simple initial conditions at 
a;o = — oo, i.e. in terms of partial differential equations with retarded boundary 
conditions. This result makes next-to-leading order computations of the average 
multiplicity feasible even in theories with strong external sources. It has possibly 
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very significant ramifications for heavy ion collisions. 

The results in section O and 0] are related in section [S] The Abramovsky- 
Gribov-Kancheli cancellations 39 , originally formulated in the context of reg- 
geon field theory models, play an important role in relating the inclusive particle 
multiplicities in ^ to moments of the probability distributions in 13 In section 
|S1 we provide a "dictionary" to map the original AGK discussion formulated in 
the language of cut reggeons into our language of cut vacuum- vacuum diagrams 
in field theories with external sources. The fact that the AGK cancellations 
require a Poissonian distribution of cut reggeons is seen to follow naturally if 
one identifies reggeons (cut and uncut) with connected vacuum-vacuum sub- 
diagrams. This identification further leads us to the conclusion that there are 
further cancellations that are not visible when the discussion takes place at the 
level of reggeons. These additional cancellations simplify the calculation of the 
multiplicity (n) as well as that of higher moments. We end this section by a 
general formula that gives the variance of the distribution of produced particles. 

In the final section, we provide a summary of the new results derived in 
this paper and of their possible ramifications for QCD at high energies. The 
paper contains four appendices. In appendix^ we give a very brief reminder 
of the Schwinger-Keldysh formalism. In appendix ^ we derive the solution 
of the Lippmann-Schwinger equation H7UI) . In appendix lUl we discuss a proof 
demonstrating that the only distributions of cut reggeons that lead to AGK 
cancellations are Poisson distributions. In appendix IdI we discuss details of the 
cancellations giving rise to the identity in eq. H109() . 

2 General framework 

In this section, we will set the groundwork for the discussion in subsequent 
sections. Much of the material here is textbook fare and can be found, for 
instance, in refs. |48l49j . However, there are subtleties that distinguish the 
discussion of field theories in the presence of external sources to discussions 
of field theories in vacuum. These are emphasized throughout, especially in 
section IT^ 

2.1 Model 

The model, for all the discussions in this paper, will be the theory of a real 
scalar field (/), with self-interactions and with the field coupled to an exter- 
nal source i{x). This simple scalar theory captures many of the key features 
of multi-particle production in QCD, while avoiding the complications arising 
from keeping track of the internal color degrees of freedom and the necessity of 
choosing a gauge in the latter. 

Without further ado, the Lagrangian density in our model is 
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The reader should note that the couphng g in this theory is dimcnsionful, with 
dimensions of the mass; the theory therefore, in 4 dimensions, is super renor- 
malizable. When we draw the Feynman diagrams for this theory, we will denote 
the source by a black dot. Keeping in mind applications of our results to the 
description of hadronic collisions in the Color Glass Condensate framework, 
one may think of j{x) as the scalar analog of the sum of two source terms 
j{x) = ji{x) + j2{x) corresponding respectively to the color currents of the two 
hadronic projectiles. 

2.2 Power counting 

Power counting is essential in understanding the perturbative expansion of am- 
plitudes in this model. For theories with weak sources, where j is parametrically 
of order 1, the leading contributions would come from the diagrams with the 
least number of vertices. We are here interested in theories (such as the CGC) 
where the sources are strong. By strong, we mean specifically that the j (f> term 
is of the same order as the other terms in the Lagrangian density. This im- 
plies that g j(x) is of order unity, as opposed to j{x) ~ 1 itself. Therefore, as 
discussed further below, independent powers of g in any given diagram arise 
only if the vertex is not connected to a source. In this power counting scheme, 
the leading diagrams can have negative powers of g. We will subsequently also 
discuss the power counting in H. 

2.2.1 Powers of the coupling 

For strong sources, gj ~ 1, attaching new sources to a given diagram does not 
change its order. There arc therefore an infinity of diagrams contributing at a 
given order in g. For instance, even the classical approximation of this model, 
which involves only tree diagrams, resums an infinite number of contributions. 
Field theories coupled to a strong source are therefore non-trivial and exhibit 
rich properties even if the coupling g is very small. 

To make this discussion more explicit, consider a generic simply connected 
diagram containing n^, external lines, internal lines, loops, n, sources 
and Uy vertices. These parameters are not independent, but are related by the 
equations, 

SUy + rij — + 2nj , 

= n, -riy -rij + 1 . (2) 

Thanks to these two relations, the order of this diagram is 

<;"vj". =5"«+2(».-l)(5j)"^ . (3) 

As explained previously, when the source is as strong as 1/g, the last factor 
is irrelevant; the order of the diagram depends only on its number of external 
legs and number of loops. Note also that the above formulas are only valid 
for a simply connected diagram: they must be applied independently to each 
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disconnected subdiagram of a more complicated diagram. Note in particular 
that a simply connected tree diagram (n^ = 0) with no external legs (n^ = 0) 
attached to strong sources is of order g~'^ . 



2.2.2 Powers of h 

It is also interesting to discuss the order in ?i of a given diagram. Recall that h 
has the dimension of an action. Therefore, 



(4) 



where Gq is the free propagator. From this equation, it is clear that the following 
mnemonic can be applied to keeping track of powers of h in the perturbative 
expansion: 

(i) a propagator brings a power of H, 

(ii) a vertex brings a power of 

(iii) a source j brings a power of l/h. 

The order in h of the generic diagram in the previous sub-section is therefore 

?i"e+"/-"v-"j = h"E+"L~'^ . (5) 

We thus recover the well-known result that the order in h, for a fixed number 
of external legs, depends only on the number of loops. For the purposes of this 
paper, it is important to note that the result does not depend on the number 
of insertions of the classical source in the diagram. Clearly this implies that 
arbitrarily complicated diagrams can be of the same order in h. In the rest of 
this paper, we do not keep track of the powers of H, and simply set h = 1. 



2.3 Vacuum— vacuum diagrams 

The forthcoming discussion in the next section, of the multiplicity distributions 
of produced particles, will be formulated in terms of "vacuum- vacuum" dia- 
grams in the presence of strong external sources. Vacuum-vacuum diagrams 
are simply Feynman diagrams without any external legs |48| . They contribute 
to the vacuum-to- vacuum transition amplitude (OoutjOin) - hence the name. 
The perturbative expansion of any amplitude with external legs generates these 
vacuum- vacuum diagrams as disconnected factors. 

When one considers a field theory in the vacuum (no external source: j = 0), 
one usually disregards vacuum-vacuum graphs because their sum is a pure phase 
that does not contribute to transition probabilities. That the sum of vacuum- 
vacuum diagrams, in this case, is a pure phase follows from the fact that no 
particles can be created if the initial state is the vacuum. The probability that 
the vacuum stays the vacuum is unity and the corresponding amplitude must 
be a pure phase. 
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The situation is very different in a theory, hke the one under consideration, 
where the fields are coupled to an external source-"^. Assume that this external 
source is such (time-dependent for instance) that an amplitude from the vacuum 
to a populated state is non-zero, namely, 



3/3^0, (/3out|0i„) ^0. (6) 



From unitarity. 



we conclude that 



2 



5] (/3out|0i„) =1, (7) 
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(Oout|0i„) < 1 . (8) 



Therefore, the amplitude (OoutjOin) cannot be a pure phase. 
It can be expressed as jJSl 

(Oout|0in) =e^^^^'l, (9) 

where iV\j\ denotes^ compactly the sum of the connected vacuum- vacuum dia- 
grams in the presence of an external source j, 

^V[J]^^Y,V^ 1^ ^k)-^ ^'^y^ ^ ■ (10) 

conn * 

In the diagrammatic expansion on the right hand side, we have represented 
only tree diagrams, but of course iV[j] contains diagrams at any loop order. 
The number preceding each diagram is its symmetry factor. Note finally that 
because external lines are absent, the power counting of the previous subsection 
tells us that vacuum-vacuum diagrams come with powers that are given by 

9'^''^-'K9jT' ■ (11) 

Therefore, in the presence of strong sources, j ~ , connected tree vacuum- 
vacuum diagrams are all of order g^^ , the 1-loop vacuum- vacuum diagrams are 
all of order 1, and so on. 

It is the squared modulus of the vacuum-vacuum amplitude that contributes 
to the transition probability. This is simply exp ( — 2ImV[j]). In the next 
section, we will discuss a method for computing 2ImV[j] directly. 

^For example, in the Color Glass Condensate (CGC) framework, incoming hadronic states 
are described by a pair of classical color sources to which the gauge fields couple. At high 
energies, such a separation is dictated by the kinematics. The resulting effective theory, the 
CGC, greatly simplifies the description of the initial state ^BI^S^S^.JS . . 

■^The factor i is purely conventional here. It has been introduced for later convenience 
when we discuss cutting rules. 
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3 Multiplicity distribution from 
vacuum-vacuum diagrams 

We will begin this section with a discussion of the cutting rules for vacuum- 
vacuum graphs in theories with external sources. These rules, for instance, can 
be employed to compute the imaginary part of the V[j] introduced in eq. (|10|) . 
We will then present an intuitive discussion of how the probability for producing 
n-particles can be expressed in terms of cut vacuum-vacuum diagrams. We 
shall demonstrate the importance of the vacuum-vacuum factors in preserving 
unitarity for a theory with strong external sources. A generating function is 
introduced to compute moments of the distribution of probabilities; we will 
show that the distribution of probabilities is non-Poissonian, even for classical 
tree level graphs. In the final sub-section, we will present a formal expression 
of the intuitive arguments developed in the previous sub-sections. 

3.1 Cutting rules 

Let us first consider the analog of the Cutkosky rules 38 in our model. These 
will also prove useful for our discussion in the following section where we develop 
an alternative method to compute multiplicity distributions. One starts by de- 
composing the Feynman (time-ordered) free propagator in two pieces according 
to the time-ordering of the two endpoints'^ 

G°(x,y) -y")G%(x,y) -xO)GO_(x,y) . (12) 

This equation defines the propagators G°_+ and G5|._. The notations for these 
objects have been chosen in analogy with the Schwinger-Keldysh formalism (see 
|5(JI51| . and the brief summary in appendix IXjl . for reasons that will become clear 
in sectional Pursuing this analogy, the Feynman propagator G° is also denoted 
as C^^. Wc introduce as well the (anti-time-ordered) propagator G^_ defined 
as 

G^__ix,y)EEeix^-y°)Gl^ix,y)+9iy"-x°)G\ix,y) . (13) 

Note that the four propagators so defined are not independent; they are related 
by the identity, 

G° + + G'L_ = G\ + G° _ . (14) 

In addition to these new propagators, we will consider two kinds of vertices, of 
type -|- or — . A vertex of type -I- is the ordinary vertex and appears with a 
factor —ig in Feynman diagrams. A vertex of type — is the opposite of a -|- 
vertex, and its Feynman rule is +ig. Likewise, for insertions of the source j, 
insertions of type -I- appear with the factor +ij{x) while insertions of type — 
appear instead with —ij(x). The motivation for introducing these additional 
vertices and sources will become clearer shortly. 

•^The superscript indicates free propagators. 
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For each Feynman diagram iV containing only "+" vertices and sources 
(denoted henceforth as iV{_| |_}) contributing to the sum of connected vacuum- 
vacuum diagrams, define a corresponding set of diagrams iV^^.y by assigning 
the symbol to the vertex i of the original diagram (and connecting a vertex of 
type e to a vertex of type e' with the propagator G^^,). Each Cj can be either of 
the "+" or "-" type. The generalized set of diagrams includes 2" such diagrams 
if the original diagram had n vertex and sources. These diagrams obey the 
so-called "largest time equation" |52I53| . If one examines these diagrams before 
the times at the vertices and sources have been integrated out, and assumes 
that the vertex or source with the largest time is numbered i, it is easy to prove 
that 

zV{...,,...} + *y{....,,..} = 0, (15) 

where the dots denote the exact same configuration of e's in both terms, for the 
vertices and sources that do not carry the largest time. This equation generalizes 
immediately to the constraint, 

where the sum is extended to the 2" possible configurations of the e's in the 
Feynman diagram of interest ^. 

This well known identity, being true for arbitrary times and positions of the 
vertices and sources, is also valid for the Fourier transforms of the corresponding 
diagrams. To see this, we first write down the explicit expressions for the Fourier 
transforms of the propagators Gj_^ 

G%{p) = - G'^_{p)= ^ . , 

pz _ _|_ jg — — le 

G'L+ip) = 2TT0{p°)5{p'' - m^) , G*[_(p) - 2tt0{-p°)5{p^ - m^) . (17) 

G^__{p) is the complex conjugate of G\^{p). Using these relations, one can 
show that 

+ - *V{+...+} + (*V{+...+})* = -2lmV . (18) 

We can therefore rewrite the imaginary part of the original Feynman diagram 
V as 

2ImF= , (19) 

where the prime in the sum indicates that the two terms where all the vertices 
and sources are of type -I- or all of type — are excluded from the sum. 

For a given term in the right hand side of the above formula, one can divide 
the diagram in several disconnected subgraphs, each containing only -I- or only 
— vertices and sources^. The -f regions and — regions of the diagram are 

*For a recent application of largest time relations to QCD, we refer the reader to Ref. 1541 . 
^Note that such disconnected subgraphs can exist only if they contain at least one external 
source. They would otherwise be forbidden by energy conservation. 
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separated by a "cut" - each diagram in the r.h.s. of eq. (|19|l is therefore a "cut 
vacuum- vacuum diagram". At tree level, the first terms generated by these 
cutting rules (applied to compute the imaginary part of the sum of connected 
vacuum-vacuum diagrams) are 

2ImV[j] = v*-f f • — • 

'^^-^^ \ ' \ \ " ^ 

+ ••• (20) 

The -f and — signs adjacent to the grey line in each diagram here indicate the 
side on which the set of -|- and — vertices is located. As one can see, there 
are cuts intercepting more than one propagator. In the following, a cut going 
through r propagators will be called a "r-particle cut" . 

It is important to note that cut vacuum-vacuum diagrams would be zero 
in the vacuum because energy cannot flow from from one side of the cut to 
the other in the absence of external legs. This constraint is removed if the 
fields are coupled to time-dependent external sources. Thus, in this case, cut 
vacuum- vacuum diagrams, and therefore the imaginary part of vacuum- vacuum 
diagrams, differ from zero. 

3.2 Probability of producing n particles 

All transition amplitudes contain the factor exp (iV[j]), which, in the squared 
amplitudes, transforms to the factor exp ( — 2ImV[j]). The power counting 
rules derived earlier tell us that the sum of all the connected vacuum-vacuum 
diagrams starts at order^ g^^. Hence, we will write 

*(Vb-]-V*[j]) = 2ImVb-]^4, (21) 

3 

where a denotes a series in g^" that starts at order n = 0; the coefficients of this 
series are functions of gj. Thus, the vacuum-to- vacuum transition probability 
is exp(— a/g^). The simplest tree diagrams entering in a/g^ were displayed in 
eq. 

We now turn to the probabilities for producing n particles. For now, our 
discussion is more intuitive than rigorous, with precise definitions to be intro- 
duced later in the section fin I3.5|l . Besides the overall factor exp(zV[j]), the 
transition amplitude from the vacuum state to a state containing one particle 
is the sum of all the Feynman diagrams with one external line. These diagrams 

^Prom now on, we shall not write explicitly the dependence in the combination gj because 
we shall assume the fields are strong and this combination is of order unity. All the coefiicients 
in the expressions we write in this section depend implicitly on gj. 
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start at order g ^. The probability to produce one particle from the vacuum 
can therefore be parameterized as 

Pi = e-«/»' ^ , (22) 

where bi is, like a, a series in 5^" that starts at n = 0. hi/g"^ can be obtained by 
performing a 1-particle cut through vacuum- vacuum diagrams. In other words, 
hi/g"^ is one of the contributions that one can find in a/g"^. Diagrammatically, 
61/3^ starts at tree level with 

+ ••• (23) 

Consider now the probability P2 for producing two particles from the vac- 
uum. There is an obvious contribution to this probability that is obtained 
simply by squaring the bi/g^ piece of the probability for producing one parti- 
cle. This term corresponds to the case wh(;re the two particles are produced 
independently from one another. But two particles can also be produced cor- 
related to each other. This "correlated" contribution to P2 must come from a 
2-particlc; c;ut through connected vacuum-vacuum diagrams. Let us represent 
this quantity as 62/5^ • Diagrammatically, 62 /^■^ is a series whose first terms are 

(24) 



The net probability, from correlated and uncorrelated production, of two parti- 
cles can therefore be represented as 



]_bi h 
2! ^ g' 



(25) 



The pref actor 1/2! in front of the first term is a symmetry factor which is 
required because the two particles in the final state are undistinguishable. 

Let us further discuss the case of three particle production before proceed- 
ing to the general case. One ("uncorrelated") term will be the cube of bi/g'^ 
(preceded by a symmetry factor 1/3!). A combination 6162/5^ will also ap- 
pear, corresponding to the case where two of the particles are produced in the 
same subdiagram, and the third is produced independently. Finally, there is 
an "intrinsic" three particle production probability corresponding to the three 
particles produced in the same diagram. We shall represent this contribution 
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by hy./g'^- More precisely, hz/g"^ is the sura of all 3-particle cuts in a/g"^. Dia- 
grammatically, some of the simplest terms in b^/g^ are 



+ ■■■ (26) 

(Only a few terms have been represented at this order, due to the large number 
of possible permutations of cuts across the various legs.) The probability of 
producing three particles from the vacuum is then given by 



1 bl 6162 b3 
3!g6+ g4 + g2 



(27) 



The previous examples can be generalized to obtain an expression for the 
production of n particles, for any n. It reads 

p=0 ^ QiH hc(p=n ^ 

In this formula, p is the number of disconnected subdiagrams producing the n 
particles, and b^/g^ denotes the contribution to the probability of the sum of all 
r-particle cuts through the connected vacuum-vacuum diagrams. This formula 
gives the probability of producing n particles to all orders in the coupling g in 
a field theory with strong external sources ^. To the best of our knowledge, 
this formula is new. Even though the dynamical details are not specified, its 
compact analytical form can be used to derive very general results for multi- 
particle production. These results will be discussed extensively in the section 
on AGK rules. Indeed, their generality follows from the fact that the dynamical 
details are not specified. A more rigorous derivation of this formula is presented 
in section EISl 



3.3 Unitarity 

In our framework, the statement of unitarity is simply that the sum of the 
probabilities P„ be unity, namely. 



E^" = i 



(29) 



ri=0 



From eq. ()28f) . it is a simple matter of algebra to show that the l.h.s. above is 
given by 



J2 Pn = e-"/''' exp 



(30) 



^Recall also that, in addition to the factors of that appear explicitly in the formula, the 
terms a, are series in g^" where, further, the coefficients of the series are functions of gj ~ 1 
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Unitarity therefore requires that 

oo 

a^J2^r ■ (31) 

This relationship between a and the br^s is in fact an identity foUowing directly 
from their respective definitions. Indeed, recall that (see eq. (|21|) ') a/g^ is equal 
to twice the imaginary part of the sum of connected vacuum- vacuum diagrams. 
Therefore, a/g^ is the sum of all the possible cuts through these connected 
vacuum- vacuum diagrams. On the other hand, K/g^ was defined as the subset 
of cut vacuum- vacuum diagrams with r-particle cuts. (Recall that these are cuts 
that intercept r propagators.) The sum of these over all values of r is therefore 
equal to a/g^ by definition. Eq. H31|) therefore confirms that the identity in 
eq. (|20|l is precisely the unitarity condition. 



3.4 Moments of the distribution 

Moments of the distribution of probabilities P„ are easily computed with the 
generating function 

oo 

F(a;)^^P„e"^ (32) 

such that (nP) — F'^p^O). Using eq. this generating function can be 

evaluated in closed form, and one obtains 

F{x) = e-°/f' exp 
The mean of the distribution of multiplicities, {n) = ^„ nPn, is 

oc 

{n)^F'{0)^-Y,rbr. (34) 

The average multiplicity is therefore given by the sum of all r-particle cuts 
through the connected vacuum- vacuum diagrams, weighted by the number r of 
particles on the cut. The second derivative of F(x) at a; = is simply 

F"(0) = (n^) . (35) 

The variance of the distribution, a = (n^) — ('^)^7 can instead be obtained 
directly from the second derivative of 

oc 

G(x) =lnF(a;) = ^y5,(e"-l) (36) 
9 t^i 



exp 



1) 



(33) 
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at X = 0. Thus, 



a = G"{Q) = \yrHr ■ (37) 
More generally, the connected part of the moment of order p reads 

("^)co„„cctcd = G(^H0)--,E'^'^^- (38) 

r=l 

One sees that if 6^ 7^ for at least one r > 1, the variance and the mean 
are not equal and the distribution is not a Poisson distribution. The converse 
is true: it is trivial to check that eq. H28(l is a Poisson distribution if foi 7^ and 
6r = for r > 1. Indeed, in this case, case, eq. H28(l would simply be 



1 fbi 



In a certain sense, cq. (|28ll is approximately Poissonian, since the term (6i/.g^)" 
has the largest power of g^^ in P„. However, even if this approximation is a 
good one for the individual probabilities P„ , it leads to an error of relative order 
unity for the value of the mean multiplicity; eq. H39|l gives = while 
the correct value is {nj = g~'^ rb^. 

3.5 Formal derivation of eq. (|28|) 

In section [3.21 the probabilities P„ for producing n particles were expressed in 
terms of br/g^. The latter are the r-particle cuts in the sum of cut connected 
vacuum-vacuum diagrams. We worked out the expressions for P„ in terms of 
br only for the cases n = 1,2,3 and guessed a generalization of the formula to 
an arbitrary n. We shall here derive more rigorously this expression for P„ and 
shall rewrite it formally to identify the br coefficients. 

An important step in the forthcoming proof is to note that the expectation 
value of the time ordered product of n fields is 

(Oo,t|r0(xi).--</>(x„)|Oin)--^...— ^e'^I^l . (40) 

This formula is equivalent to the usual way of obtaining time-ordered Green's 
functions from a generating functional. However, in our case, the source j is not 
a fictitious variable that one sets to zero at the end but is instead the (physical) 
external source. We can also use the standard LSZ reduction formula 03 to 
rewrite the amplitude for producing n particles as * 



(Pl---P„out|Oin) = ^ j 



i5j{xi) 



e'^I^'l . (41) 



*The wave- function renormalization factors Z correspond to self-energy corrections on the 
cut propagators of the vacuum-vacuum diagrams. 
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The probability for producing n-particles is given by the expression 



n (2^)32^;, 



(42) 



where Ei = \J ip^ + rn? . Substituting the r.h.s. of eq. in the above, and 
noting that 



is the Fourier transform of the expression in eq. (|17fl , we can write the probabihty 
Pn directly as 



P. 



tV[j+] ^-tV*[j-] 



3+=3~ 



where j_] is the operator 



Z 



6 



Sj+{x) 5j-{y) 



(44) 



(45) 



The sources in the amplitude and the complex conjugate amplitude are labeled 
as j+ and j_ respectively to ensure that the functional derivatives act only on 
one of the two factors. 

In eq. H44|l . the mass m in G\_ and in □ + is the physical pole mass of 
the particles under consideration rather than their bare mass. The role of the 
prefactor Z~^ becomes more transparent if we write 



V ■ 



d'^xd'^y ZGl_(x,y) 



a, 



5j+{x) Sj-{y) ' 



(46) 



and recall that Z is the residue of the pole of the renormalized propagator. All 
the factors under the integral in eq. (|46|l are therefore renormalized propagators, 
or their inverse. 

Using eq. (|44|l . we can write the generating function F(x), defined earlier in 
eq. (|22I), as 

F{x) = e'^'^IJ+'^-l e'^[J+I e-'^*[J-l . (47) 

j+=j-=3 

The reader may note that the factor exp(iV[j+]) represents the sum of the 
usual Feynman vacuum- vacuum diagrams with source j+. Likewise, the fac- 
tor exp(— iV*[j_]) is comprised of the complex conjugate of vacuum-vacuum 
diagrams, with source j_. One thus begins to see a formal realization of the 
intuitive discussion of the previous section in terms of cut diagrams. Alter- 
nately, the second factor of eq. (|47|l can be interpreted as containing terms in 
the Schwinger-Keldysh expansion 50 51] that have only -I- vertices, while the 
third factor contains terms that have only — vertices. A brief introduction to 
the Schwinger-Keldysh formalism is provided in appendix 1X1 
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We shall now further elaborate on the interpretation of P„ (and F{x)) as cut 
vacuum-vacuum Feynman diagrams as well as their relation to the Schwinger- 
Keldysh framework. Consider the action of the operator on the 

vacuum-vacuum factors in eq. (|T7|l . The operator (O^ + m'^)6 /Sj+{x) takes 
a diagram of exp(iV[j+]), removes one of its sources, and amputates the prop- 
agator to which this source was attached. (This procedure includes any self- 
energy decoration on this propagator because of the factor Z~^.) Z~^(\I\y + 
m?)5 / 5j-{y) does the same on the factor exp(— iV*[j_]). Finally, the factor 
ZG^_{x,y) is a (dressed) propagator that links the two diagrams that have 
been opened by the functional derivatives with respect to the sources. When 
a; = 0, one can convince oneself that the action of the operator exp(I?[j-|_, j_]) 
is to build all the other vacuum-vacuum diagrams of the Schwinger-Keldysh 
formalism, namely those that have at least one G'^_^ or G'1_, 

Here iVc;i^[j+, j-] denotes the sum of all connected vacuum- vacuum diagrams 
in the Schwinger-Keldysh formalism, with the source j+ on the upper branch 
of the contour and likewise, j- on the lower branch. When the magnitudes of 
the sources on the upper and lower branches are the same, j+ = j- = j, it is 
well known that this sum of all vacuum- vacuum diagrams is equal to 1. (This 
statement is proved in appendix 1X1 ) Therefore, eq. H47() satisfies F{0) = 1, as 
unitarity dictates it should. 

When X ^ 0, the generating function F{x) is the sum of all vacuum-vacuum 
diagrams in the Schwinger-Keldysh formalism, with the two "off-diagonal" prop- 
agators substituted by 

G^. + ^ e^G% , 

Gl_ -> e^G°._ , (49) 

and where the source is the same on the upper and lower branches of the contour. 
Naturally, the logarithm of this generating function is obtained by keeping only 
the connected vacuum- vacuum diagrams in the sum. 

Recall now that the building blocks of the Schwinger-Keldysh perturbative 
expansion are the same as those obtained by the cutting rules for calculating 
the imaginary parts of Feynman diagrams. We can therefore identify each term 
in the expansion of F{x) with a cut vacuum-vacuum diagram. Moreover, one 
sees from eq. ()49|l that the power of in a given term is equal to the number 
of propagators G'l_+ or G° _ it contains. Therefore, the order in equals the 
number of particles on the cut. Since \n{F{x)) contains only the cut connected 
vacuum-vacuum diagrams, and since br/g^ has been defined earlier as the sum 
of all the r-particle cuts through connected vacuum- vacuum diagrams, we have 
just proved that 

oo , 

ln(F(x)) = -4+^e-^. (50) 
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The coefficients br can thus be formally identified and computed in principle. 
One obtains the constant term in this formula, — a/g^, from the limit x — > —oo 
of eq. I@7I) and from eq. (EH • 



4 Direct calculation of the average multiplicity 

In the previous section, we obtained a general expression for a generating func- 
tion which can be used, in principle, to compute all moments of the multiplicity 
distribution. In this section, we will concentrate on the first moment, the aver- 
age multiplicity. We will develop the Schwinger-Keldysh formalism, introduced 
previously in section 13.51 and in appendix ^ to compute this quantity both 
at leading and next-to-leading order. While the result at leading order is well 
known, the next-to-leading order result is new, and is the central result of this 
paper. 



4.1 General formula 

Although one could in principle compute the average multiplicity using eq. (|34|) , 
it would be extremely impractical because this method requires one to evaluate 
separately all the r-particle cuts through the vacuum- vacuum diagrams^. It 
is much more straightforward to obtain the multiplicity from the generating 
function given in eq. H47|) . We obtain 



3+=3-=] 



where we have used eq. (^SJ- Writing out exphcitly using eq. 

this expression can be rewritten as 



(51) 

3, 



(n)^ J d^xd^y ZGl^{x,v) — 



(52) 



3+=0- 



Functional derivatives of the sum of connected vacuum-vacuum diagrams are 
connected Green's functions with the number of external points equalling the 
number of derivatives. When we differentiate iVo^, 



5 



(53) 

the corresponding Green's functions are the path-ordered Green's functions de- 
fined in eq. H113I) of appendix^ More precisely, they are the connected com- 
ponents of these Green's functions - hence the subscript "conn" . The operators 



^This would be an almost impossible task even at tree level. 
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Z~-^{\I\ + m?) in eq. (|52ll amputate the external legs of these Green's functions 
(along with the self-energy corrections they may carry). 

Denoting r'^'^i "^")(a;i, • • • , Xn) as the amputated version of the Green's func- 
tion in eq. (|53|l . we finally express the average multiplicity as^" 



(n) 



d^xcfy ZG%_{x,y) 



rW(x)r(-)(2/)-i-r(+-)(x,y)' 



(56) 



This formula is valid to all orders in the coupling constant and is a key result 
of this paper as will become clearer later. Diagrammatically, we can represent 



it as 



{n) 



(57) 



where the shaded blobs represent the amputated Green's functions T. In prin- 
ciple, one should amputate as well the self-energy corrections attached to their 
external points and use the full — h propagator. Alternatively, one can use 
the bare — h propagator, and keep the self-energy corrections attached to the 
external legs of the shaded blobs. 

Since the F's are connected Green's functions, the second of the two diagrams 
is at least a 1-loop diagram. Therefore only the first diagram can contribute at 
leading order in the coupling. The functions T are obtained with the standard 
Schwinger-Keldysh rules, described in the appendix^ Note also that, unlike the 
expressions for the probabilities, there are no vacuum-vacuum diagrams left in 
this formula. They have all disappeared thanks to the fact that exp(iVs^ [j, j]) = 
1. This makes the evaluation of the multiplicity much simpler than that of the 
probabilities P„ even for n = 1. 



4.2 Leading order result 

The multiplicity at leading order, namely 0(g~^(gj)"), is obtained from the left 
diagram in eq. (|57|l with each of the two blobs evaluated at tree level. We have. 



^"if one is only interested in the calculation of the multiplicity (n), a more elementary 
derivation is to start from 

<«) = / (Oin|aLt(p)aout(p)|Oin) . (54) 

A reduction formula for the correlator that appears under the integral gives 

(Oin|aL(p)<*out(p)|Oin) = ly d^xd'^ye-'^-^e'^-y 

y.{U:, + rr?){Uy + 'rr?) (0in|9i(a;)<^(j/)|0in) . (55) 

The expectation value on the right hand side, in wrhich the fields are not time-ordered, is just 

the propagator G y{x,y) of the Schwinger-Keldysh formalism. Splitting it into connected 

and disconnected contributions, one can see that eqs. lEU-OSJ and jSU arc equivalent. 
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where the sum is over all the tree diagrams on the left and on the right of the 
propagator _ (represented in boldface) as well as a sum over the labels +/— 
of the vertices whose type is not written explicitly. At this order, the mass in 
G'^_ is simply the bare mass, and Z = 1. 

Eq. (|58|l can be expressed very simply in terms of the retarded solutions of 
the classical equations of motion. The identities, 

Lr + + - + _ - G-^ , 

G*L+ -G"__=Gl , (59) 

where is the free retarded propagator, will be essential in demonstrating this 
result. Consider for instance, one of the tree subdiagrams contained in eq. I|58|l . 
Summing over the + or — indices of the "leaves" (or "blobs") of the tree, the 
outer layer of propagators is transformed into retarded propagators, thanks to 
eqs. (jSnj- These then give factors of / d'^yGii{- ■ ■ , y)j{y) that are independent 
of the indices of the vertices immediately below. This trick is applicable here 
because the source j{x) is identical on both branches of the closed time path. 
By repeating this procedure recursively until one reaches the "root" of the tree, 
one turns all the propagators into retarded propagators. As a consequence, one 
finally obtains the same tree diagram (ending at point x or y), where all the 
propagators are now retarded propagators. 

Now recall the perturbative expansion of the retarded solution of the classical 
equation of motion^ ^ 

iD + m^)Mx) + ^4{x)^j{x), (60) 

with an initial condition at = ~oo such that the field and its derivatives 
vanish 

lim (j)cix) = , lim d^'cjjdx) = . (61) 

x*^ — > — oo — * — oo 

In that case as well, the solution can be expressed, identically, recursively, in 
terms of the retarded Green's function. 

The sum over all leaves up to the root of the tree on each side of the cut in 
eq. can therefore be identified as 

(62) 




Note that in eq. (|58|l . we need to amputate these tree diagrams, so it is (□ + 
m'^)(j)c{x) that enters in eq. H58|) rather than (pdx) itself. More precisely, we 
have 

(n) f d''xd^yGl_{x,y) (n^ + m^){ay+m^)cf>,{x)My) ■ (63) 



-'-'^To see that the various factors of i present in the Feynman rules of the Schwinger-Keldysh 
formalism agree with this identification, simply multiply eq. l6Ui by i, and recall that we have 
defined the propagators as the inverse of «(□ + m^) - for instance, see eq. 1171 . 
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Using the explicit momentum space form for G'J__ in eq. 143|l .we can rewrite 

2 



(n) 



(64) 



This leading order result is well known. In the CGC framework for gluon produc- 
tion, the multiplicity was evaluated by computing the spatial Fourier transform 
of the classical field 0c at late times p6i.27.30.31.33.34] . To connect their work 
to our eq. we use the identity 



\d^o + EDMx) = do (e'P-^ [do ~ lEp] M^)) 



(65) 



and perform an integration by parts of the r.h.s. of eq. Il64|) (using the initial 
condition of eq. iGl\i ) to obtain 



s'P-^ (□ + m^)(j)cix) = lim / d^x e'?"^ [d^„ - iEp] ^x) 



(66) 



In practice, the r.h.s. of this equation is the preferred method for evaluating the 
average multiplicity. One only needs to solve the classical equation of motion 
with the initial conditions of eq. I)61|l and perform a spatial Fourier transform 
of the solution at the latest time. 



4.3 Next-to-leading order result 

At next-to-leading order (NLO) in the coupling constant {0{g'^ {gj)")), one ob- 
tains two sorts of corrections to the leading contribution of eq. H58|) : 

(i) the right diagram in eq. (|57|l with the blob evaluated at tree level 




(ii) 1-loop corrections to diagrams of the kind displayed in eq. (|58|l . These 
come from the left diagram of eq. (|57|l . with one of the blobs evaluated at 
1 loop, and the other blob kept at tree level 




In both cases, these contributions provide corrections of order g'^ to the multi- 
plicity (n), to be compared to the leading contributions of order g^^. 



20 



That one obtains these two sets of contributions at NLO is also weh known. 
However, what is novel, and will be demonstrated below, is that both these 
contributions can be computed entirely with retarded boundary conditions and 
retarded propagators. This makes it feasible to perform NLO computations 
(to all orders in the external sources) by solving equations of motion for the 
classical and small fluctuation fields with boundary conditions entirely specified 
at t = —oo. 

4.3.1 Evaluation of eq. (|67|l 

We begin by discussing the diagrams displayed in eq. (|67|l . These topologies 
are well known since they are the ones involved for the production of fermion- 
antifermion pairs 55 . (Indeed, for this process, they constitute the lowest order 
contribution.) Evaluating the diagrams of eq. (|67|l is equivalent to calculating 

the tree level propagator attached to an arbitrary number of tree diagrams 

of the type depicted in eq. 1)62(1 . Each of these attachments, thanks to eq. ((62|l . 
can be replaced by the classical field (pc itself. Thus 



where a dotted line terminating in a cross represents an insertion of the classical 
field (eq. H62|l ). At the vertices where (pcS are inserted, the +/— indices must 
still be summed over; a vertex of type — simply corresponds to a sign change for 

the field insertion. This is equivalent to computing the propagator G-i (a;, y) in 

the presence of the background field (j)c. We shall mimic the method used in |55j 

and write down a Lippmann-Schwinger equation whose solution is {x,y)- 

For later reference, we shall write it for any component Gee'ix, y) as 



G,Ax,y) = GlA^.y) -^9 G%ix,z) </>c(^X,- G,,Az,y) , (70) 



where the resummed propagator has no superscript 0. cr^ = diag(l, —1) is the 
third Pauli matrix, and its purpose in eq. (|7Uf) is to provide the correct minus 
sign when the vertex at which the classical field is inserted is of type — . 

The Lippman-Schwinger equation is solved by performing a simple rotation^^ 
of the +/— indices. The details of this procedure can be followed in appendixIBl 
The solution of the Lippmann-Schwinger equation for the "rotated" propagator 
G reads 



Such transformations were introduced in )56) and discussed more systematically in )57l . 




(69) 




(71) 
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where the resummed retarded and advanced propagators are given by 

oo oo 

G. - Gl [<^K] " > = ] " , (72) 

n=0 n=0 

and 

G,^G,{GI)-'gI{GI)-'g,. (73) 

Here and are the free retarded and advanced propagators respectively, 
and G° = G% + G°__ = 2tt5{p^ - m^). The combination G^G^ is the full 
retarded propagator in the presence of the background field (j)c, amputated of 
its rightmost leg. The propagators G^^' of the Schwinger-Keldysh formalism are 
obtained by performing the inverse of the rotations described in appendix IbI and 
the explicit expressions are given there. 

To compute the next-to-leading order contribution to the multiplicity (n) 
from eq. ijHTji . we need specifically the propagator G+_. The solution of the 
Lippman-Schwinger equation gives 

G+- = G,Gl-'G'l_Gl-'G, (74) 

We are not done yet because eq. (|56f) tells us that one has to amputate the 

external legs of G-| This requirement is conveniently fulfilled by introducing 

"retarded and advanced scattering T-matrices" , defined by the relations 

R — + ^R^R^R ' 

G, ^G° -HG°r,G° . (75) 
The contribution of eq. (|67|) to (n) is then given by 

jj^0^J d'xd'ye-^^-e^^y j d\d'vT,{x,u)G\_{u-v)TAv.y) • 

(76) 

In momentum space, it reads more compactly as 

^">-o =/ (4%/ (4%^"-^'''"'^^^^"''''^ • ^^^^ 
Using the identity that r^(— g,p) — [T^(p, — g)]*, we obtain finally, 

This formula is the equivalent for scalar bosons of eq. (50) in and eq. (94) 
in [HHI- 

The retarded T-matrix Tji{p, —q) can be computed from the retarded solu- 
tion of the partial differential equation , 

{U + ra^ + g(l)c{x))T^{x) , (79) 

l3\Ye remind the reader that <; is a dimensionful quantity. 
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where 77 is a small perturbation about the classical background field (l)c', the 
equation is therefore the linearized equation of motion in this background. Using 
Green's theorem^'' (see appendix A of for an extended discussion), 



r){x) = « / d yG^{x,y) dy„ rj{y) , 

Jj^o=const 



as well as eq. (I73|) . one can prove that 



lim 



[9^0 - iEp] 7](x) 



(80) 



(81) 



where r]{x) is the retarded solution of eq. (|79|l with initial condition ri{x) — e**'^ 
when xq ~oo. One should not confuse this equation with cq. despite 
the fact that they look similar. When using eq. I|()t)|) one must solve the classical 
equation of motion in the presence of the source j in order to obtain In 
eq. (|81|l . one is instead solving the equation of motion for a perturbation 77 
about the classical solution c/fc- The eigenvalues of this perturbation arc used 
to construct the small fluctuations propagator in the background field. The 
analog of eq. (|81|) for fermion production was previously used in j58| and |36l37j 
to evaluate the yield of quark-antiquark pairs in pA collisions and in the initial 
stages of nucleus-nucleus collisions respectively. 



4.3.2 Evaluation of eq. (|68|l 

We will now consider NLO contributions to the multiplicity from 1-loop cor- 
rections such as those depicted in eq. (|t)8|) . Only the subdiagram that contains 
a closed loop is new. Using the same trick as in eq. (|66f) . the contribution of 
eq. I)68ll to (n) can be written as 



(n) 



(2) 

NLO 



d^p 



{2'KY2Ep 

X 



lim 

Xq — * + 00 

lim 

Xq— *-t-oo 



d^x e''P-^[9o -«-Bp]0e(a;) 
d'x e'J'-''[do-iEp]d)cAx) 



(82) 



The novel quantity considered here is the 1-loop correction to (l)c{x) - we will 
represent this quantity as (pcii^)- We will discuss below our strategy to compute 
this quantity by solving partial differential equations with boundary conditions 
at x^ —00. 

All the diagrams contributing to 0c, 1 (a;) include a closed loop, to which may 
be attached an arbitrary number of tree diagrams. All but one of these tree 
diagrams are terminated by sources j; this last one has an external leg at the 

^■'The prefactor i is due to the fact that the inverse of □ + is iG according to our 
conventions for the propagators. 
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point a;, 



+/- 




(83) 



The first simplification one can perform is to use eq. H62|l to collapse all the 
trees terminated by sources into insertions of the classical field 0c. This gives 



+/- 



( 

X 




9 


+ 

i 


( 



o 



(84) 



Again, at each remaining vertex or insertion of 0c, one must sum over the type 
+/— . The number of insertions of 0c can be arbitrary, both on the loop and 
on the leg leading to the final point x. Therefore, eq. H84|) can be written more 
explicitly as 



(85) 



where all the propagators in this formula must be evaluated in the presence of 
the background field (pc- 
We can rewrite 0c. i as 



-«5 



/ d^y G++{y,y)^tG+^{x,y) 



ig I d'^y G^{x,y) G++{y,y) , 



(87) 



where we have used eq. (|59f) . This expression is equivalent to the retarded 
solution of the following partial differential equation: 

(□ + m^ +50c(a;))0c,i(a;) = -gG++{x,x) , 



^^From eqs. 11231 of appendix IbI we see that the propagators G++{y,y) and G iy,y), 

which enter in eq. I85i . are both equal to 



G++iy,y) = G—{y,y) = -G ^G^-^G^G^^-^G 



(86) 
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with an initial condition such that 0c, i and its derivatives vanish at = — oo. 
The source term G'++(x,x) in eq. is given by eq. 



d'^ud'^v 
d'^q 



0-1 



{x,u)G"{u,v) 



G"r^G, 



2TTS{q^ - m^) / d-^u 



'0-1 



d% 



G"r^G, 



{v,x) 
ix,u) e*«-" 
{v,x)e~^^-\ (89) 



using 

From the r.h.s. of the relation 



(2^ 



2nSiq^ - n?) . 



(90) 



d^u 



0-1 



(x, u) e*«'' = lim i / 0^(2.^^)9^^ e*'?" = ?7(x), (91) 



we see that this quantity is precisely the retarded solution at point x of eq. H79|l 
with an initial condition ri{u) = e"? " at uq = — oo. Further, by noticing that 



{v,x) 



GrGI- 



{x,v) , 



(92) 



the factor involving G^~^G in eq. H89|) is obtained in a similar way, except that 
the initial condition must be e~*^ ". 

Following this discussion, we can now outline here an algorithm for obtaining 
the quantity (l)c^i{x) needed in the calculation of 

(i) For an arbitrary on-shell momentum g, find the solutions •q'q^^x) and 



r}q [x) of eq. (|79|) whose initial conditions at a;o ^ are respectively 
ri['^\x) = e*9-^ and r^['\x) = e"^?-^ 

(ii) Calculate 6*++ (a;, x) as 

d^q 



G++{x,x) = \ J A^2nSiq'-m')^l+\x)vi-\x) 
(iii) Solve eq. H88() with a vanishing initial condition at = — oo 



(93) 



For completeness, we note that G'++(a;,x) contains ultraviolet divergences. 
They appear as a divergence in the integration over the momentum q in eq. (|93() . 
By going back to the diagrams that are contained in G++(x, x), it is easy to see 
that these divergences are the usual 1-loop ultraviolet divergences of the (j)^ field 
theory in the vacuum, namely, the tadpole and the self-energy divergences^^. 

^'^Of course, the scalar toy model with a <f>^ self-interaction that we are considering in this 
paper is somewhat peculiar in the ultraviolet sector, since this theory is super-renormalizable 
in four dimensions. But the procedure outlined here for removing the divergences would be 
valid in QCD as well. 
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These divergences correspond to the two terms of eq. in which the loop is 
dressed by zero or one insertions of the classical field (pc (the terms with two 
insertions of 4>c or more on the loop involve the 1-loop corrections to the n-point 
functions for n > 3, which are finite in four dimensions). From this observation, 
it is clear that one must subtract the following quantity from G++(x,a;) 

SG++{x,x)=5i + {5^n + S„,2)(l),{x) , (94) 

where Si is the counterterm that makes the tadpole finite, and and the 
two counterterms that make the self-energy finite. One possible approach is to 
evaluate both eq. 193|) and the counterterms with the same momentum cutoff 
A, and to subtract eq. (|94|l from eq. H93|) before letting A go to infinity. 

5 AGK cancellations 

We will discuss here a set of combinatorial rules that relate the discussions 
of the n-particle probabilities in section |21 and the average multiplicity in the 
previous section. We find that these combinatorial relations are identical to 
the Abramovsky-Gribov-Kancheli (AGK) cancellations first discussed in 
the context of reggeon field theory. This enables us to map out explicitly, in 
the context of multi-particle production, the terminology of field theories with 
external sources (such as the Color Glass Condensate) with that of reggeon 
field theory. Our discussion is in fact more general than the original AGK 
discussion, since it is not limited to diagrams that have an interpretation in 
terms of reggeons alone. 

Before we proceed, we shall provide here a very brief sketch of the extensive 
literature on the subject of the AGK cutting rules. A review of early work 
can be found in ref. 001 ■ Applications of AGK rules to reggeon field theory 
models of particle production in proton-nucleus and nucleus-nucleus collisions 
were discussed at length by Koplik and Mueller More recently, several 

authors considered the application of these rules in perturbative QCD. For a 
comprehensive recent discussion, see Ref. [421431 and references therein. We 
note that AGK rules have also been considered previously in the Mueller dipole 
picture 01] of saturation/CGC |45I46| . A detailed discussion and additional 
references can be found in |T7| . 

5.1 AGK-like identities 

The AGK relations are a set of identities that relate the average multiplicity 
of particles produced in hadronic collisions to the multiplicity of a single cut 
reggeon. Thanks to these identities, all the other contributions to the average 
multiplicity involving more than one reggeon exchange (where the reggeons that 
are not cut correspond to absorptive corrections) cancel in the calculation of the 
average multiplicity. 
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Even though our discussion is based on a completely different theory, one can 
obtain a set of identities that play the same role as the original AGK identities^'^. 
To see these, consider Pn,m, the term that has a factor in the expression 

eq. (EH for P„, 

^-"-i E E (95) 

p+q=m ai-\ \-ap=n 

Obviously one has 

oo 

Pn - ^"'"^ ■ (96) 

Here m is the number of disconnected subdiagrams contributing to P„; each 
disconnected diagram starts at order g~^. In eq. H95() . (7 is the number of dis- 
connected subdiagrams that are not cut - they do not contribute to the number 
n of produced particles, and correspond to the "absorptive corrections" of |39| . 
Likewise, p is the number of subdiagrams that are cut - they contribute to the 
number of produced particles. In the model of |39j. m would be the number of 
reggeons in the diagram. 

There is however one important difference between our considerations and 
those of ISni- In our paper, n is the number of produced particles, while in 
n is the number of cut reggeons, regardless of how many particles are produced 
in each of them. Therefore, in our case, n can be larger than m (because each 
cut subdiagram can produce many particles), while n is always smaller or equal 
to m in 

We shall now introduce the term of order l/g'^"\ Fm{x), in the generating 
fimction F{x) introduced in section 3. It is defined as 

00 

Pm{x^ — ^ ^ Pn,ni ^ 

This formula can be obtained either by using eq. H95I) and performing the sum 
explicitly, or by expanding eq. to the appropriate order 1/5^™. Evaluating 
the successive derivatives of Fm{x) at a; = 0, we can easily prove the identities 

if m > 2 , 
if m > 3 , 



To prove these formulas, we used the unitarity relation, a — br- The physical 
interpretation of these formulas is quite transparent. They tell us that Feynman 

^^In this sense, one can understand these identities to be a very general consequence of field 
theories with strong external sources. 



/ 00 \ ™ 

J2 ^e"" - a . (97) 



nPn,Tn = , 

n=l 

00 

n{n - l)Pn,ni = , 



(98) 
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diagrams that have two or more disconnected subdiagrams do not affect the 
mean multiplicity, that Feynman diagrams with three or more disconnected 
subdiagrams do not affect the variance, and so on... These identities are the 
generalization to our situation of the eqs. (24) of 

Note that the sum over n must be extended to infinity due to the above 
mentioned fact that we are counting particles instead of cut reggeons. 

5.2 Closer connection with reggeon theory 

One can in fact make more detailed connections with the discussion of ref . 139) , 
and in a sense, generalize the results there. Recall that the latter is formulated 
in the language of reggeons so the connections to our framework are not obvious. 
We remind the reader that in eq. (|28|l . for P„, the index p represents the number 
of cut subdiagrams contributing to the production of the n particles. Therefore, 
one can further divide P„ and write the probability for having an event in which 
n particles are produced in p cut subdiagrams^'* 

PiCl-^^'^^ E (99) 

QiH [-ctp=n ^ 

which of course is related to P„ by 

n 

Pn=Y.Pni- (100) 

The unspecified "subdiagrams" in our discussion are a generalization of the 
"reggeons" of Since the discussion in ,39, was at the level of the cut 

reggeons rather than at the level of the individual particles, we may obtain 
a closer correspondence to the results in |39| by integrating out the index n 
representing the number of produced particles. Summing over n in eq. H99|l 
gives the probability of having p cut subdiagrams, 

n—p ^ Vi* / 

Clearly, the multiplicity of cut subdiagrams has a Poissonian distribution. This 
is not a surprise because they are disconnected from one another^^. Moreover, 
the average number of cut subdiagrams is equal to 

^*This new quantity Pn'^p should not be confused with Pn,m- The latter is the probability for 
producing n particles in m subdiagrams, irrespective of whether cut or uncut. Defining Pn,p,q 
as the probability of producing n particles with p cut subdiagrams and q uncut subdiagrams, 
the former quantities are respectively P„,m = Y.p+q=m Pri,p.q and p}{^l = Y.^^ Pn,p.q- 

^^In this result was achieved thanks to an assumption about the impact factors, while 
in our model it follows naturally because particle production arises only through the coupling 
to an external source. 



28 



One may define the average number of particles produced in diagrams with p 
cut subdiagrams as^" 

(,) =pgr^. (103) 



n.p 



This resuh implies that 



{n)^p{n)^. (104) 



In other words, a diagram with p cut subdiagrams produces on average p times 
the number of particles produced by a diagram with a single cut subdiagram. 
Further, one can deduce immediately that the average number of produced par- 
ticles, (n) = <7~^ X^r average number of cut subdiagrams multiplied 
by the average number of particles produced in one cut subdiagram, 

(n) = (neut) (n)^ . (105) 

This property was assumed in ref . (HHI ; they therefore focused only on properties 
of the distribution of cut reggeons. In our approach, this property arises natu- 
rally, regardless of the nature of the subdiagrams we consider. (This result is 
therefore more general than the ladder diagrams that correspond to reggeons.) 

In order to exactly reproduce the identities in [^21, we introduce the proba- 
bility TZp^m of having p cut subdiagrams for a total number m of subdiagrams 
(m — p of which are therefore uncut). Expanding the exponential in eq. ()101|l 
to order m — p, one obtains 

1 1 / a 

{m-p)\pl V g'^J Vf^ 



This distribution of probabilities then satisfies the relations 

rn 

if m > 2 , ^ p'R-p,m = , 



p=i 



if m > 3 , ^Pip- l)7?.p,rn = , 

p=2 



(107) 



This set of identities is strictly equivalent to the eqs. (24) of [23 • For instance, 
the first relation means that diagrams with two or more subdiagrams cancel 
in the calculation of the multiplicity. These relations are therefore a straight- 
forward consequence of the fact that the distribution of the numbers of cut 
subdiagrams is a Poisson distribution. They do not depend at all on whether 

^''This quantity must be defined as a conditional probability. Indeed, the denominator is 
necessary in this definition in order to take into account the fact that events with different 
numbers of cut subdiagrams can occur with different probabilities. 
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these subdiagrams are "reggeons" or not. In appendix [CI we prove the reverse 
result: the only distribution TZp for which cqs. (|107|l hold is a Poisson distribu- 
tion. 

It is also instructive to note that these identities obeyed by the distribution 
of cut / uncut disconnected subdiagrams have been obtained by "integrating out" 
the number n of produced particles. By doing so, one goes from the distribution 
H99|l - which still depends on all the quantities 61, 627 ^3 ■ ■ • ~ to H101() . in which 
all the dynamics is summarized in one number a = br- Integrating out the 
number of particles has the virtue of hiding most of the details of the field theory 
under consideration. This is in fact the reason why the AGK identities have a 
range of validity which is much wider than the reggeons models for which they 
had been derived originally. 

However, this robustness has a cost: by integrating out the number of parti- 
cles, one has lost a lot of information about the details of the theory, which means 
that certain questions can no longer be addressed. For instance, in eq. I|1Q4|I , the 
factor {n)^ - the average number of particles produced in one cut subdiagram - 
can only be calculated by going back to the microscopic dynamics of the theory 
under consideration. 



5.3 Additional cancellations in ^n) 

In the original approach of Abramovsky, Gribov and Kancheli, one would esti- 
mate the particle multiplicity predicted by the reggeon model as follows: 

(i) compute the diagram with one cut reggeon exchange in order to obtain 
the average number of cut reggeons, (ricut) 

(ii) calculate the average number of particles in one cut reggeon (?^)-^. This 
second step in general requires a more microscopic description of what one 
means by "reggeons" 

The advantage of our approach is of course that it is entirely formulated as a 
field theory in which the elementary objects are particles. This makes it possible 
to compute these two quantities from first principles. Indeed, expressions for 
these in terms of cuts through vacuum-vacuum Feynman diagrams are 

('^cut) = , {n) = ^ . 108 

As one can see, in our microscopic description, there is a further cancellation in 
their product, between (ricut) and the denominator of (n)^. This result could 
be anticipated because we already know that the average multiplicity is equal 
to g~'^Y.r'^br. 

Moreover, by comparing eqs. H34|l and 1)58(1 . we see that we have the following 
identity at leading order, 
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In the r.h.s. of this relation, the types of the endpoints of the propagator repre- 
sented in boldface must be held fixed, and one must sum over the topologies and 
over the +/— indices for the trees on both sides of this propagator. The dia- 
grammatic expansion of the left hand side of eq. (|109|l is of course obtained from 
eqs. (1231, (ESJ, if^ - etc... We also know from eq. that the tree diagrams 
that appear in the r.h.s. are nothing but the retarded solution of the classical 
equation of motion. In other words, eq. 1] 109(1 tells us that the sum of all cuts 
through (time-ordered) tree connected vacuum-vacuum diagrams, weighted by 
the particle multiplicity on the cut, can be cast into a much simpler expression 
in terms of retarded tree diagrams. This means that additional cancellations 
must occur in the left hand side of eq. (|109|l in order to have this dramatic 
simplification. These cancellations are discussed explicitly in appendix IdI 



5.4 Higher moments 

The cancellations discussed in the previous subsection, which are in fact respon- 
sible for the identity of eq. H1U9|) . are crucial in order to be able to calculate 
the particle multiplicity in a simple way. We have seen explicitly that one can 
calculate the multiplicity at leading order directly from the classical field 0c ) 
and in principle at higher orders as well. Of course, the NLO calculation has 
the complications of a typical 1-loop calculation; in particular one has to treat 
the ultraviolet divergences. 

It turns out that higher moments of the multiplicity distribution are also 
calculable from 0^, by formulas that are tractable even though they are more 
complicated than in the case of the multiplicity. Let us consider the case of the 
variance as an illustration. It is obtained as the second derivative at a: = of 
the log of the generating function, lii{F{x)). Therefore, it is given by 

a ^ (n^) - {nf = [PbV, J-] + 2?'^+, J-]] e^""^- | • (HO) 

connected 

Performing the functional derivatives and dropping the disconnected terms that 
show up gives 



(j = {n)+ J d'^xd'^yd'^ud'^v G\_{x,y)G\_{u,v) \T^+-+-\x,y,u,v) 

+V'^-++\y,u,v)T^+\x) + T^++-\x,u,v)T^^\y) + {x ^ u,y ^ v) 
-f r(++) [x, u)T'-—\y, v) + r(+-) (x, i;)r(-+) (y, u) 

+r(+-) (x, 'y)r(-) (2/)r(+) (u) + r(-+) (y, M)r(+) (a;)r(-) (w) 

-hr(++)(a;,w)r(-)(y)r(-)(w) +r(— )(y,v)r(+)(a;)r(+)(w)j , (111) 
where the F's are amputated Green's functions in the Schwinger-Keldysh for- 
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malism. Diagrammatically, the variance can be expressed as 



(112) 

All the diagrams represented in this equation are deviations from a Poisson 
distribution (for which we would simply have a = (?^))- The diagrams of the 
first line start at tree level, while those of the second and third lines start 
respectively at 1 and 2 loops. Similar formulas can be derived for even higher 
moments, but they become increasingly complicated as the order of the moment 
increases, even at tree level. 

6 Summary and Outlook 

In this paper, we studied particle production in a model scalar field theory 
coupled to a strong (j ~ l/g) time dependent external source. This theory is a 
toy model of the Color Glass Condensate formalism for particle production in 
high energy hadronic collisions in QCD. We shall summarize here some of the 
novel results of this paper. 

• We obtained a general formula (eq. (|28|l 'l for the probability P„ to produce 
n particles in a field theory with time dependent external sources. (A 
formal expression of this formula is given in eqs. (|44|l and (|45|) .'l 

• We derived the corresponding generation function for moments of the mul- 
tipUcity distribution. We showed that the generating function (eq. l|T7|) ') 
can be related to the sum of all vacuum-vacuum diagrams in the Schwin- 
ger-Keldysh formalism with off-diagonal propagators given by eq. (|49|l . 

• We obtained a general formula (eq. (|56|l ) for the average particle multi- 
plicity in terms of the amputated Schwinger-Keldysh Green's functions. 
While the result to leading order is well known, the result at next-to- 
leading order (NLO) is new and remarkable. It is not surprising of course 
that the average multiplicity at NLO can be expressed in terms of small 
fluctuation fields. What is non-trivial and remarkable is that these fields 
can be computed by solving small fluctuation equations of motion with re- 
tarded propagators and initial conditions specified at x° — > — oo (eqs. ()79f) . 
iP|l and P^). 

• We outlined an algorithm to compute these NLO contributions. Such an 
algorithm has been implemented previously for fcrmion pair production in 
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the Color Glass Condensate, where the diagram of cq. Ijt)?!) is the leading 
contribution. It has not been implemented for the diagram of eq. (|68|l . 

• We observed that the Abramovsky-Gribov-Kancheli cancellations, origi- 
nally formulated in terms of relations among reggeon diagrams, can be 
mapped into the field theory language of cut vacuum-vacuum diagrams. 
Disconnected vacuum- vacuum subgraphs can be identified as reggeons and 
corresponding cut sub-graphs as cut reggeons. 

• We showed that an immediate consequence of this mapping is that AGK 
cancellations in reggeon field theory hold if and only if the distribution of 
cut reggeons is Poissonian. This follows because two disconnected cut sub- 
diagrams do not interact by definition! The conjecture is proved explicitly 
in appendix 10 We observed that there are additional cancellations (as 
represented by eq. (|109|l ) that arc not visible at the level of the AGK 
cancellations. 

• We derived a general formula for the second moment (the variance) of the 
multiplicity distribution. This formula can be expressed in terms of the 
average multiplicity plus "non-Poissonian" terms. These are nonzero even 
at tree level in a field theory with strong time dependent external sources. 

We shall now discuss sonic of the ramifications of our results, in particular 
for high energy QCD. Our studies were indeed motivated by their possible rel- 
evance to hadronic collisions at very high energies. In particular, in the Color 
Glass Condensate approach, hadronic collisions at very high energies are de- 
scribed by an effective field theory where the degrees of freedom are wee parton 
fields coupled to strong time dependent external valence sources. A key feature 
of the CGC approach is that weak coupling techniques are applicable. The tech- 
niques developed here can therefore be applied to study particle production^'^. 
Of particular interest is particle production in the earliest stages of heavy ion 
collisions, and the possibility that one can understand the thermalization of 
a quark gluon plasma from first principles. As is immediately clear from the 
formalism developed here, the solution of such problems is intrinsically non- 
perturbative (even in weak coupling) requiring the summation of infinite series 
of Feynman graphs. Numerical algorithms performing these summations are 
therefore essential. 

Indeed, first steps in carrying out the program outlined here have already 
been performed. Gluon production 29 30 31 32 33 34 and quark pair produc- 
tion |36I37| have been computed numerically previously to leading order in the 
coupling and to all orders in the source density. There are strong hints that a 
deeper understanding of thermalization requires an analysis incorporating en- 
ergy loss effects [6 1162163] . These may include the scattering contributions of the 
bottom-up scenario |64I65| , or equivalently the energy loss induced by collective 
effects |66I67I68] . In the CGC framework, these effects can be shown to appear 

However, one should keep in mind the fact that transition amplitudes in a slowly varying 
background might be difficult to define rigorously )60| . 
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at next-to-leading order in the coupling. Qualitative studies |()9I35| of the NLO 
effects suggest that these may be very important in driving the system towards 
equilibrium. Our study here is a first step towards a quantitative formulation 
of this problem. In a follow up study [70], we will address the connections of 
our approach to kinetic approaches |71I72| . 
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A Schwinger-Keldysh formalism 
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Figure 1: The closed time path used in the Schwinger-Keldysh formalism. 

Our intuitive discussion of cutting rules for computing probabilities in field 
theories coupled to external sources, begun in section 12.31 and continuing in 
section 131 culminates in the formal expression of these, in section 13.51 in terms 
of the Schwinger-Keldysh vacuum-vacuum amplitudes. The Schwinger-Keldysh 
formalism was developed in the 1960's in the context of many-body quantum 
field theory |50I51I| . It is used primarily in thermal field theory problems in order 
to calculate thermal averages of operators [Z3|- Since it may be less familiar to 
those working outside this area, we will provide a brief review of the formalism 
before proceeding to prove a result stated in section 1^51 

The primary feature of the Schwinger-Keldysh formalism, relative to the 
standard formulation of quantum field theory, is that the time fiow of paths is 
defined on a contour comprised of two branches (see the figure^ wrapping the 
real axis. The upper branch, denoted +, runs along the real axis in the positive 
direction while the lower branch, denoted — , runs in the opposite direction. The 
aim is to compute the correlation functions 

G^^--">(X1 • • • a;„) = (Oin|Pc c/^^^'H^l) ' ' ' 4>^'''H^n)\0in) , (113) 

where Pc denotes a path ordering of operators along the time contour. Operators 
are ordered from left to right by decreasing values of their curvilinear abscissa 
on the contour, with operators living on the lower branch placed to the left of 
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operators living on the upper branch. When restricted to the upper branch, 
it is equivalent to the usual time ordering. Similarly, it is equivalent to an 
anti-time ordering on the lower branch. In eq. (|113|l . the indices can have 
the values +/— indicating the branch of the contour where the corresponding 
field is located. This definition of the path ordering suggests that the correlator 
defined in eq. H113|) can also be written as 

G{^--">(xi---x„) = (Oi„|T[ J] cf>^'^\x,)\T[ II 0(^^)(a;,)]|Oi„) . (114) 

Note that the vacuum state used in order to define this correlation function is 
the "in" vacuum state on both sides^^. 

The perturbative calculation of the correlator in eq. (|113|l requires special 
Feynman rules. These are obtained in the Schwinger-Keldysh formalism and 
their derivation is analogous to that of regular time-ordered correlation func- 
tions. One first expresses the Heisenberg fields in terms of the free field 0in of 
the interaction picture. One can then rewrite eq. H113() as 

G^^— >(a;i • • •a;„) = {0,^\Pc <j)t\^i) ■ ■ ■ ^t\xn) g/c '''-^.-(^-C-)) |0i„) . 

(115) 

Note that the time integration in the exponential now runs over the full con- 
tour C. The perturbative expansion of this correlator is obtained by expanding 
the exponential. This procedure differs from the usual Feynman rules in the 
following two ways: 

• at each vertex, the time integration runs over the full contour C 

• the free propagator connecting vertices is the path-ordered product of two 
free fields 

G\x,y) = (Oi„|Pc'^i„(a;)</'i„(y)|Oi„) . (116) 

These rules are implemented by first breaking the time integration at each vertex 
in two pieces. One corresponds to an integration over times in the upper branch 
while the other integrates over the lower branch. The vertices come in two 
varieties, -I- vertices and — vertices, and the indices must be summed over 
at each vertex. For example, the contribution of the — vertices is written as 
an integration over the real axis in the positive direction modulo the addition 
of a minus sign to the vertex of type — . This follows from the fact that the 
integration runs in the negative direction on the lower branch of C Likewise, the 
propagator in eq. H116|l is broken into four different propagators, corresponding 
to where the two fields lie on the contour C, 

Gl+ix,y) = (Oi„|T(/)i„(a;)0i„(y)|Oin) , 
G°_(a;,y) = (Oi„|T0i„(a;)</>i„(2/)|Oin) , 

G°+(x,2/) EE (Oin|0in(a;)0in(y)|Oin) , 

G°_(X,2/) = (Oin|<^in(y)</'i„(x)|Oi„) . (117) 



^^This feature is a consequence of the fact that the thermal average of an operator is a 
trace of the operator muhiphed by the density matrix. 
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Note that the Fourier transforms of these free propagators are precisely those 
encountered previous in the derivation of cutting rules (see eq. (fTTj) '). The 
propagator G°^, connects a vertex of type e to a vertex of type e'. 

In the discussion in section 13.11 the ± vertices, propagators and sources 
were introduced as a formal trick to compute the sum of the Feynman vacuum- 
vacuum amplitude and its complex conjugate amplitude. We now understand 
these as arising naturally from the Schwinger-Keldysh formalism which captures 
the physics of theories in external fields/finite temperatures. For example, an 
identity (eq. that we showed in section ITTl is equivalent to a property of 

the Schwinger-Keldysh formalism used in section [3.51 the sum of all vacuum- 
vacuum diagrams is equal to one. It remains true even if the field is coupled 
to an external source, provided that source is the same on the upper and lower 
branches of the contour C. 

The proof of this result in the Schwinger-Keldysh formalism is very simple. 
Consider a connected vacuum- vacuum diagram, and integrate out all but one of 
its vertices. The result is a function of a single point we call x. This function, by 
virtue of not possessing external legs, and being coupled to a source that has the 
same value on the -I- and — branches of the contour, has the same magnitude on 
both branches of the contour. Therefore, when we sum over the +/— type of this 
last vertex, we get an exact cancellation because the — index has the opposite 
sign and the same magnitude as -I- index. Thus all connected vacuum-vacuum 
diagrams are zero when the sources are identical on both legs of the contour. 
Because the sum of all vacuum- vacuum diagrams, is the exponential of the sum 
of the connected ones, it is therefore equal to unity. 



B Solution of Lippman-Schwinger equation 

The Lippman-Schwinger equation H70|) can be solved by performing the following 
rotations of the propagators and field insertions: 

Gee' Gap = UaeUpe'Gee' 

e,e'=± 

Kf^^ J2 UaeUpe'cri, (118) 



£,e'=± 

with 

Under these rotations, the free matrix propagator and field insertion matrix 
become 

and G° are the free retarded and advanced propagators, and G^ stands 
for the combination 

Gl = G%+G°__. (121) 
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Its Fourier transform is G'°(p) = 2tt5{p'^ — m?). In terms of these new prop- 
agators and vertices, the Lippmann-Schwinger equation has the same form as 
that of eq. H7U|) . modulo the replacements of eq. (|118|l . The main simplification 
follows from observing that the product GS'^ is the sum of a diagonal matrix 
and a nilpotent matrix. This makes it easy to calculate its n-th power. In 
particular, we have^^ 

( \ 
G°[e3gT= ■ (122) 

Note that the off-diagonal equations mean that in the new basis the Lippmann- 
Schwinger equations do not mix the retarded and advanced propagators. 

The result above can be easily inverted to obtain the propagators of the 
Schwinger-Keldysh formalism. One obtains^"* 

y^lO — 1 y^O y^O — 1 

+ - — '-'R^R 

G++ = - [Gj^G" ^G"G° ^G^ + (G„ -I- G^)] 

G- = \ {G,Gl-^GlG\-^G, (G„ + G J] . (123) 

As one can see, all these propagators resumming the effect of the classical back- 
ground field 0c can be expressed in terms of the retarded and advanced prop- 
agator in the background field. This is a useful property because retarded 
propagators can be obtained by solving a simple partial differential equation 
with retarded boundary conditions. (Advanced propagators can be obtained 
from the retarded ones by permuting the endpoints.) 



C Reciprocal of eqs. ( 11071 ) 

It is also instructive to see what the AGK identities, as given in eqs. H1Q7|I . 
imply for the distribution of multiplicities of cut subdiagrams in eq. H1U1|) . Let 
us assume the following generic form for the probability 7?.p of having p cut 
subdiagrams : 

Tl^ = R(T^)-T^Tf , (124) 

^•^This formula is true even if each propagator has different arguments. Therefore, it is also 
true for the convolution of n -|- 1 propagators, between which one inserts the classical field. 
■^*One can check that this solution satisfies the following properties 

G++ -HG__ = G_+ + G+_ 

G++ - G+_ = , G++ - G_+ = G^ . 
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where 77 is a parameter that counts the number of subdiagrams (in the situation 
studied in section lS?^ the quantity a/g"^ played this role). The prefactor R{ri) is 
determined from the requirement that the sum of these probabilities be unity, 



R{il) = 



+00 ^ 

lp=Q 



(125) 



and the coefficients rp must be positive in order to have positive probabilities. 
Moreover, we can assume without loss of generality that ro = 1. The H-p^q that 
appear in eqs. (|107|l are obtained by expanding the prefactor R{ri) to order q—p, 
and are thus completely determined from the coefficients rp. 

The relations p07|l are very complicated non-linear relations when expressed 
in terms of the coefficients rp, because of the prefactor R{rf). In order to inves- 
tigate their consequences, it is much easier to introduce the function 

+00 _! 

G(x)^^-rp(77x)P. (126) 
p=0 ^' 

Note that the prefactor R{ri) is nothing but G~^(l). Therefore, if we consider 
the following generating function for the probabilities TZp, 

+00 

H{x) = ^npxP , (127) 

p=0 



it is related to G{x) via : 
Recalling now that 



H{x) = G{x)/G{l) . (128) 



^P=II^P,9. (129) 

q=p 

and using eqs. H1U7|) . we have the following relations for the successive derivatives 
of H{x) at X = 1 : 

i^(l)(l) = 7^l,l , i^(2)(l)^2!7^2,2, ••• h^^\\) ^ p\np^p . (iso) 

Therefore, we can rewrite the generating function H{x) as 

+OC +00 
p=0 p=0 ^ 

In order to obtain the second equality, we have used the fact that Rp^p is obtained 
by keeping only the order in the expansion of the prefactor R{ri), which is 
equal to unity from our choice to have vq = 1. Therefore, the AGK identities, 
eqs. (|107|) . imply that 

Vx, G{x) = G{x - 1)G{1) . (132) 



+OC +00 

H{x) = J2 ^P.p(^ - 1)'' = II - - 1)" = G{x - 1) . (131) 



38 



One sees immediately that the exponentials, G{x) — exp(aa:), are trivial 
solutions of this functional relation. Note that this solution corresponds to 
a Poisson distribution for the probabilities TZp. We are now going to argue 
that this obvious solution is in fact the only solution which is compatible with 
the positivity of the coefhcients rp. Let us factor out of G{x) the exponential 
behavior it may contain by writing 

G{x)=g{x)e^\ (133) 

where g{x) is a function whose growth at infinity is bounded by a polynomial. 
The functional equation (|132|l . when rewritten in terms of (/(x), implies 

Vx, ff(x + l)=g(x)g(l), (134) 

and in particular 

VneZ, g{n) = g{0)g'\l) . (135) 

Since g{x) cannot grow faster than a polynomial at infinity, we must therefore 
have g{l) = 1, and cq. H134|l becomes 

Vx, g{x+l)^g{x). (136) 

Hence the function g(x) must be a periodic function of period unity. This is 
where the positivity of the rp's plays a role: it implies that all the successive 
derivatives of G{x) are positive for x > 0. In terms of g{x), the n-th derivative 
of G{x) reads 

n 

G(") (x) = e"^ 9^"'^ (x) , (137) 

p=0 

where the = n\/p\{n — p)\ are the binomial coefficients. We must therefore 
find a periodic function ^(x) such that 

n 

Vn e N , Vx , ^ CP a"-'' ^(^^(x) > .. (138) 

Let us write the periodic function g{x) as a Fourier series, 

<?(x)^5]ff,e2-'^ (139) 

where g-i = g* since ^(x) is real when x is real. The left hand side in the 
inequality (|137|l can be written as 

5oa" + 2^|g,| {a^ + Air^f)''/^ cos{2nlx + n(j)i) , (140) 
1=1 

where we parameterize gi{a+2iTTl) = \gi \ vt? + 47r^Pexp («</);) . We see that for 
any n and I, we can choose a value of x such that the cosine equals —1. Moreover, 
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by going to sufficiently large n, we can make (a^ + 47r^Z^)"/^ arbitrarily large 
compared to a". Therefore, the values of \gi\ such that eq. 1] 140(1 is always 
positive have to be chosen arbitrarily small. Therefore, we have proved that the 
only choice of the Fourier coefficients which is compatible with the positivity of 
the probabilities TZp is 

yiel* , 91^0 . (141) 

Thus, the function .9(2;) is constant and equal to 1 (from g(l) ~ 1), and the 
only functions G{x) that obey the functional relation H132|l are the exponentials 
G{x) = exp(aa;). We have therefore proved that the only distributions TZp that 
lead to the AGK cancellations are Poisson distributions. 



D Discussion of the cancellations in eq. ( 11091) 

We will here study the cancellations in cq. H109|l more explicitly. Note first that 
the terms in 17"^ rhr must correspond to the terms in the r.h.s. of eq. ((109|l 
topology by topology. For the simplest topology that appears in hi/g"^, only 
1-particle cuts are allowed, and it is trivial to check the identity of eq. H109() . 
The first non-trivial topology is the "star" diagram that appears in bi / and 
in h^l g^ ■ Its contribution to the left hand side of eq. H109|l is 




On the other hand, the contribution of this topology to the right hand side of 
eq. (|109|l is comprised of the terms 

Note that we have discarded two vanishing terms, in which a vertex of type + 
(resp. — ) was completely surrounded by sources of type — (resp. +), because 
such a configuration is forbidden by energy configuration (because of the 0(±po) 

factors in the Fourier transform of the propagators G |_ and G^ ). At this 

point, it is a trivial matter of inspection to check that eqs. ((142|l and H143|l are 
identical. One should stress that it was crucial to weight the various r-particle 
cuts hrjg^ by the multiplicity r on the cut. 

One can also see that the equivalence stated in eq. ((109|) between the method 
of calculating (n) from the classical field and the direct method implies that 
62 ^ (in other words, this equivalence fails if one disregards the 2-particle 
cuts in eq. (|142|l 'l. Incidentally, this proves that the distribution of produced 
particles is not a Poisson distribution, according to the discussion at the end of 
section 13 
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